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Detuning-controlled internal oscillations in an exciton-polariton condensate 
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We theoretically analyze exciton-photon oscillatory dynamics within a homogenous polariton gas 
in presence of energy detuning between the cavity and quantum well modes. Whereas pure Rabi 
oscillations consist of the particle exchange between the photon and excitons states in the polariton 
system without any oscillations of the phases of the two sub-condensates, we demonstrate that any 
non-zero detuning results in oscillations of the relative phase of the photon and exciton macroscopic 
wave functions. Different initial conditions reveal a variety of behaviors of the relative phase between 
the two condensates, and a crossover from Rabi-like to Josephson-like oscillations is predicted. 

PACS numbers: 71.36.-l-c, 67.85.Fg 


Exciton-polaritons are new mixed eigenmodes result¬ 
ing from strong coupling between the photon state in a 
microcavity and exciton state in a quantum well, which 
inherit properties of both light and matter [l[. Polari- 
tons interact without dephasing due to elastic collisions 
in their excitonic component, while the photon compo¬ 
nent provides them with extremely light effective mass 
allowing Bose-Einstein condensation (BEC) at high crit¬ 


ical temperatures^ 
of experimental 



. In the recent ye ars , gre at wealth 

works in 


and theoretical [E 
the field of polariton BEC demonstrated the coherent ef¬ 
fects analogous to those in atomic condensates 


superconductors [H, 341, or liquid Helium At 

the same time, due to their non-trivial dispersion, short 
lifetime, and therefore non-equilibrium nature of con¬ 
densation, polariton quantum properties differ con¬ 
siderably from those appearing in BECs of atoms or in 
superconducting systems. The phenomena investigated 
include condensation in traps [j, H, superfluidity 
El , yortices ssEia solitons (l3. [S polariton 


poiari 

Wk 


olarization features |2ll-l23l|. and Josephs on p henomena 


24l-l28j . Latest experiments to date [39|, l40| report 


obseryation and hne control of such inherent features of 
polariton systems as relaxation oscillations and Rabi os¬ 
cillations. There has as well been some theoretical effort 
dedicated to polariton Rabi dynamics, reyealing possible 
ways to increase the coherence time and proposing qubits 
and logic gates based on exciton-polaritons H 4^. and 
stochastic processes within polariton condensates j43l] . 


Strong coupling regime which is considered here is de¬ 
fined as the Rabi splitting between the npper (UP) and 
lower (LP) polariton branches at the anticrossing being 
large compared to their emission linewidths [l| . The fre¬ 
quency of this two-leyel oscillator depends on the Rabi 
frequency UIr and the detuning between energies of bare 
photon and exciton states at zero wayevector. From the 
point of yiew of the initial exciton and photon states, 
if the detuning is zero, the polariton gas is a half-and- 


half mixture of the two coherent components constantly 
performing mutual transformations, hence one may say 
that Rabi oscillations are density oscillations between the 
‘photon’ and ‘exciton’ condensates. The picture becomes 
more complicated when the photon and exciton disper¬ 
sions are shifted with respect to each other: internal oscil¬ 
lations between exciton and photon components involye 
oscillations of the relative phase of the two macroscopic 
wave functions. 

In this work, we present fully analytical investigation of 
internal oscillations in the two-component polariton sys¬ 
tem and discuss possible regimes of the dynamics. We 
consider an idealized polariton gas with constant chemi¬ 
cal potential and neglect all non-equilibrium effects asso¬ 
ciated with particles gain and dissipation. For this model 
condensate, we demonstrate that at any non-zero detun¬ 
ing different types of oscillations are possible, from har¬ 
monic and anharmonic modifications of Rabi oscillations 
up to the transition to a so-called Josephson regime anal¬ 
ogous to internal Josephson effect in a two-state BEC of 
^’^Rb atoms 3, ilj . We also address interaction-induced 
corrections to our analytical solutions brought into the 
system by possible increase of the polariton density. 

Within the mean field approach [l^ . temporal evolu¬ 
tion of the macroscopic wave functions of cavity photons 
ifc and quantum well excitons 'ifx is described by the 
coupled Schrbdinger and Gross-Pitaevskii equations. 
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with Eq ^ the bottoms of energy dispersions and mc,x 
the effective masses of photons and excitons. g > 0 is the 
constant of repulsive exciton-exciton interaction. Parti¬ 
cle transfer between the subsystems is described by the 
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coupling term ^ hnji/2, and we neglected any external 
potentials and spin degree of freedom. 

When the system is in the strong coupling regime, the 
polariton state is an eigenstate with an equal (in the ab¬ 
sence of interactions) superposition of a photon and an 
exciton. The positive sign chosen in Eqs. ([T|), ([2]) in front 
of the coupling term hflii/2 imposes that the antisym¬ 
metric mode ('i/'c — 4’x)/y/^ with the relative phase tt 
is the lower energy level {i.e. corresponds to the LP 
state) while the symmetric mode {ipc + '>px)/V2 with 
zero relative phase is the upper one. An initial state of 
the polariton system, being some linear combination of 
these two modes, results in density oscillations between 
the photon and exciton subsystems. When interactions 
are present, the effective lower energy level is blueshifted 
(while the upper energy level appears redshifted), and 
the eigenmodes are no longer the antisymmetric and the 
symmetric ones. Still, the relative phase oscillations dis¬ 
cussed below go around the time-average values tt or 0. 

Considering homogeneous case for simplicity and as¬ 
suming momentum equal to zero, we neglect spatial 
derivatives in Eqs. © and ©■ After the transforma¬ 
tion V'C.x {t) = \/nc,x{t) we get four non-linear 

dynamical equations for photon and exciton populations 
nc,x{t) and phases Sc,x{t): 

dtnc,x = Ty/ncnx sin [Sc — Sx) , (3) 

cos (S'c - S'x), (4) 

2\j nc 

dtSx = -e°x - gn-x - cos {Sc - Sx), (5) 

2\j nx 

where we have rescaled lengths and energies in terms of 
yJfijmc^R and respectively, time as t^R —?► t, and 
the wave functions as ipc,x /\/^ /mc^R '4’C,x- 

In order to investigate the dynamics, it is conve¬ 
nient to introduce new variables: relative phase S{t) = 
Sc{t) — Sx{t) and population imbalance p{t) = {nc{t) — 
'nx{t))/n, where n = nc{t)+nx{t) is the total number of 
polaritons. Variables p and S obey the coupled equations 

p = - \/\- p^ sin S, (6) 

• gngn p 

The dimensionless detuning <5 = Cq — Cx blueshift 

value gn/2 are the parameters which determine different 
regimes of the system behavior. Eor a system with con¬ 
stant chemical potential the equations ©, © are Hamil¬ 
tonian with the conserved energy iJ(5', p) = {6—gn/2)p+ 
gnp^/^+yJl - cos 5, where total population n is con¬ 
stant. Equations © and © admit analytical solution in 


terms of quadratures: 


cos S = 


TT _ gn 
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After obtaining formal solution of the evolution equa¬ 
tions, it is worth noting that the interaction constant g is 
of the order of 10“^ (estimated from g = 0.015 meV-/im^ 
i) , while the unsealed n is of the order of unity (which 
corresponds to ~ 10^° cm“^ Q). Hence, for the closed 
conservative system the blueshift parameter gn/2 is al¬ 
ways of the order of 10“^, and the effect of interactions 
on internal oscillations should be negligible. We esti¬ 
mate the upper limit for realistic polariton densities as 
Ucr ^ lO^^-lO^^ cm“^ [46| . In our numerical simulations, 
we raised the total density n up to 0.5 x 10^^ cm“^. How¬ 
ever, as shown further, even for large densities neglecting 
the interactions results in little loss of accuracy. 

In the absence of interactions, the integral in © re¬ 
veals an explicit solution for the population imbalance: 


P{t) = 


hS 


Va;2 - 


sin(a;t — (p). 


( 10 ) 


where h = (5p(0) -I- — p(0)2 cosS'(O) = const is the en¬ 

ergy (per particle) defined by the detuning and the initial 
conditions, w = is the renormalized frequency of 

internal oscillations (in scaled units, it corresponds to 
+ 5^), and (p = arcsin[(/ii5 — p(0)w^)/-\/a;^ — h^]. 
The relative phase between the photon and exciton sub¬ 
systems is then given by © with the substitution of (HOD . 

The phase-plane portrait of the conjugate variables 
{p, S) is shown in Eig. ©a)-(c) for three detuning val¬ 
ues. When the detuning compensates the blueshift {6 = 
gn/2)^ for any initial conditions the system performs fi¬ 
nite motion along the selected trajectory (depending on 
the energy). This case is displayed in Eig ©a) and it cor¬ 
responds to Rabi-like oscillations. Without interactions, 
if the initial state is prepared in such a way that the ex¬ 
citon and photon populations are equal, i.e. p{0) = 0, 
the system stays in the pure Rabi regime of density oscil¬ 
lations with the time-average (p) remaining zero in time 
and without any change of the relative phase S = tt. Any 
non-zero p(0), however, will result in harmonic oscilla¬ 
tions in both population imbalance and relative phase 
with the Rabi frequency Hr {lo = 1). Allowing for in¬ 
teractions in this case results in the frequency given by 
a; = 1 -I- gnh/2 [d^. The obtained correction reveals the 
decrease of frequency with the amplitude {h < 0), which 
is characteristic for nonlinear pendulums, however in this 
case the frequency does not reduce to zero at the separa- 
trix (/i = 0), tending instead to its non-perturbed value 
Hr. 
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FIG. 1. (Color online) (a)-(c) Phase-plane portraits of the 
conjugate variables p and S for different detunings. 
Trajectories of different colors correspond to different values 
of the energy h according to the scale given on the right, (a) 
5 = (?u/2, (b) 5 = —0.5, (c) 5 = —1.5 (all energies are in the 
units of hQii). For positive detunings, portraits flip with 
respect to the vertical axis. The red dashed and black 
dotted lines represent separatrices for particle densities 
10^° cm~^ and 0.5 x 10^^ cm“^, respectively, (d)-(f) 
Relative phase S = Sc — Sx against time, for the 
trajectories of the phase-plane portrait (b) as marked. Right 
to left: (f) anharmonic modification of Rabi oscillations, (e) 
sawtooth-like oscillations at the separatrix, (d) the regime of 
running relative phase (internal Josephson effect). For 
gn/2 = 0.002, analytical expressions (1101) . ([8| coincide with 
direct numerical solutions of Eqs. ([51), (O (the solid lines). 
The dashed lines show the numerical solutions for 
(;n/2 = 0.5 x 10“^. Other physical parameters: 
me = 0.6 X 10“^mo, hflji = 5 meV, g = 0.015 meV-pm^. 


Fig- mb) andme) show the phase space trajectories for 
small and large negative detunings, respectively. In this 
case, the two-component system can evolve in two differ¬ 
ent dynamical regimes, depending on the initial energy. 
Closed trajectories representing finite motion at low en¬ 
ergies belong to the regime of Rabi-like oscillations sim¬ 
ilar to the previous case. The difference consists of the 
renormalization of the oscillation frequency w, the anhar- 
monicity of the relative phase oscillations (see Fig.mf))i 
and the shift of the time-average population imbalance 
to a non-zero value h5/uP'. This regime of oscillations 
is kind of an interplay between the modified Rabi dy¬ 
namics and an analog of internal Josephson effect: for 
small-amplitude oscillations, one may say that the shift 
of natural frequency corresponds to Josephson “plasma 
frequency” LOjp = Sflp. 

As the energy h increases at fixed detuning (or, alter¬ 
natively, as I (51 increases at fixed h), the phase oscillations 
grow in amplitude up to 7r/2 and acquire strongly anhar¬ 
monic sawtooth profile shown in Fig. me), while the tra¬ 
jectory on the phase-plane portrait approaches the sep¬ 
aratrix line defined by cos S' = 6y /(1 — p)/(l + p)- After 
crossing the separatrix, one sees a dramatic change from 
Rabi-like to Josephson-like dynamics: while the density 


imbalance oscillates around its new equilibrium value 
/iJ/vTT^, the relative phase between the photon and 
exciton condensates S(t) becomes monotonously increas¬ 
ing (or decreasing, depending on the sign of S) in time 
like shown in Fig. mb)- This regime of running phase is 
analogous to the a.c. Josephson effect in the Josephson 
junction 36|, or to internal Josephson-like oscillations be¬ 


tween the populations in a mixture of spin-up and spin- 
down atoms when external magnetic field is applied j44l | . 
(N.B., all the above explanations imply that the system 
is that of lower polaritons. For upper polaritons, the 
relative phase would oscillate around zero instead of tt, 
and the decrease instead of increase in h would bring the 
system closer to the separatrix and consequently to the 
Josephson regime). 

Numerical solutions of Eqs. © and © taking into 
account interactions start to noticeably differ from the 
analytical solutions (HOD and © owhen the dimension¬ 
less parameter gn/2 becomes of the order of 10“^ and 
larger (see Fig. mb)-(f))- First-order analytical correc¬ 
tion to the separatrix lines appears as cos S' = (J — (1 — 
p)gn/4:)y/{l — p)/{l -I- p), and is shown in Fig. m8')“(c) 
for gn/2 = 0.05 as black dotted lines. The region of fi¬ 
nite motion corresponding to modified Rabi oscillations 
is slightly reduced by interactions, although for higher 
detunings their influence weakens. 

Let us address the preparation of initial states, namely 
p(0) and S(0), which define the energy h. When LP 
and UP branches are resonantly excited by two spectrally 
narrow, phase-correlated laser pulses, the effective state 
created in the system is a linear combination of the lower 
and upper polariton states with controllable populations 
and relative phase. For example, if the pulses arrive in 
phase, the initial state will be purely photonic, and if 
they arrive in antiphase, it will be purely excitonic. The 
multitude of intermediate cases produce all the variety of 
possible initial conditions for the considered system. 

All the dynamical regimes are finally summarized in 
the energy-detuning diagram displayed in Fig. [2] In the 
absence of interactions, dimensionless energy of the LP 
system (energy per particle) can change in the range 
—Cl) < h < 0, while for upper polaritons 0 < h < uj. 
Corrected by interactions. 


h - in in iH --n ^ 

^min/max r^ + 2 (jj' 2 


1 ,,,, 
2):? + —I- <“> 


The critical values of h which correspond to the tran¬ 
sition between the ‘modified Rabi’ and the ‘internal 
Josephson’ regimes are defined for each detuning as h = 
±((5 — gn/2) + gn/A. The diagram h{S) is divided in four 
regions which correspond to Rabi-like and Josephson-like 
oscillations of lower and upper polaritons, as shown in 
Fig. [2] The regime of pure Rabi oscillations (with con¬ 
stant relative phase) corresponds to the point S = gn/2, 
h = —1. (Or, for the hypothetical equilibrium UP con¬ 
densate, 6 = gn/2, h = -1-1.) It can also be seen from 
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FIG. 2. (Color online) Energy-detuning diagram 
summarizing the internal dynamics (in the units of hnn). 
Minimal and maximal energies (per particle) are given by 
the solid lines ~ tVI + as marked. For LP, 

h < 0 and (S) = tt, and for UP, h > 0 and (S) = 0. Within 
these limits, the different dynamical regimes are divided by 
the separatrices h = ±5. The dashed lines take into account 
interactions for the polariton density 0.5 x 10^^ cm“^ and 
g — 0.015 ‘R’-regions correspond to modified 

Rabi oscillations. ‘J’-regions correspond to Josephson-like 
oscillations with the running relative phase. Red dotted line 
is a guide to the eye for the corresponding data in Fig. [S] 


this diagram that the larger is the detuning between the 
modes, the less extra energy is needed for the transition 
to the Josephson regime to happen. 

Finally, we analyze (for vanishing interactions) how the 
described internal photon-exciton dynamics influence the 
phase of the photon field. Using the solutions for p{t) and 
S{t) given by (ITUl) and ([5]), we analytically integrate Eqs. 
O and ©• The solutions read: 


Sc,xit) = 5'c,x(0)- ^ t 


— arctan 


— arctan 


(w2 ± hS) tan ± - /i2 

uj {h± S) 

(tu^ ± hS) tan ^ =F y/uj'^ — h'^ 


to {h± S) 


■ ( 12 ) 


This result highlights the fact that linear rotation of the 
photon and exciton phases given by — (^C + £'y)l/2 is 
modulated by additional periodic terms. Fig. [3] shows 
the periodic parts of the phases (na for different initial 
states at fixed detuning One can clearly see 

that when the polariton system is in the regions ‘R’ or 
‘J’ (as marked in Fig. la and far from the separatrices, 
the additional periodic terms in both Sc and Sx are 
approximately harmonic with frequency w and of small 
amplitude (see Fig. [SKa) and (i)). Upon approaching the 
separatrix h = —S, the additional oscillations of the pho¬ 
ton phase become strongly anharmonic (Fig. |3Kb)~(e)), 


and their amplitude grows up to tt/2. Fourier spectrum 
of the periodic part of Sc then consists of multiple fre¬ 
quencies. Being added to the linear term, these sawtooth 
oscillations result in the ladder-like behavior of the pho¬ 
ton phase. At the same time the evolution of the exciton 
phase Sx stays practically unchanged. On the contrary, 
when the energy is close to another critical value h = -|-5, 
the exciton phase experiences the ladder-like behavior 
while the photon phase oscillations are close to harmonic- 
ity (see Fig.[3Ke)-(h)). Experimentally, the photon phase 
and its evolution can be determined through interferom¬ 
etry with an external reference with a well-defined 
phase. Even though the photon phase alone is not enough 
to retrieve the whole information about the relative phase 
5, it can clearly identify the crossover between the two 
different dynamical regimes while the initial conditions 
are being changed. 

It is important to note that, when taking into account 
the out-of-equilibrium nature of polaritons such as their 
dissipation and replenishing of the system from the reser¬ 
voir, one expects that the system should show relaxation 
oscillations to the stable points shown in FiglT](a)-(c). 
However, this appears to be the case only only when the 
gain and loss rates are linear and constant. Modelling the 
dynamics with non-linearities such as gain saturation 271 
or the reservoir dynamics results in new regimes of 
evolution, with the detuning and the initial conditions no 
longer playing a crucial role, therefore they are subject 
to a separate study. 

In conclusion, we analytically analyzed the influence 
of the photon-exciton energy detuning and repulsive in¬ 
teractions on the internal oscillatory dynamics of the po- 
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FIG. 3. (Color online) Periodic parts of the photon and 
exciton phases Sc,xit) as marked, for S — 0.5 and different 
initial conditions according to the red dotted line in Fig. [2] 
(a) h = 0.95/ijjjjjj, (b) h = 0.6h ^;„ , (c) h = 0.448/ijjjjjj, 
(d) h = 0.44hjjjjjj, (e) h = 0, (f) h = 0.44/imax, 

(g) h = 0.448/imax, (h) h = O.Ghmax, (i) h = 0.95/imax. 
The background color identifies the regimes (‘R’ or ‘J’) in 
agreement with Fig. 0 For negative detunings the curves 
corresponding to the photon and exciton phases swap. 
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lariton system. We demonstrated that at any non-zero 
detuning, the two-component system can, depending on 
its energy, oscillate around its equilibrium point or tran¬ 
sit to the regime of monotonously growing relative phase, 
which we connect to internal Josephon effect. Interac¬ 
tions, on the contrary, are shown not play a qualita¬ 
tive role in the presented physics. While the Joseph- 
son regime we describe is very much analogous to that 
of conventional bosonic Josephson junctions, the signif¬ 
icant difference is that the effects we report lie in the 
non-interacting regime of the exciton-polariton system. 
When present, interactions do not significantly modify 
the dynamics, in contrast to the situation in strong¬ 
interacting atomic systems, where interactions could un¬ 
veil the regime of macroscopic quantum self-trapping of 
populations [ 2 ^. At last, we predict that the crossover 
between the two regimes of dynamics can be experimen¬ 
tally observed by detecting the photon phase evolution 
in photoluminescence from the cavity. 
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ful discussions. This work is partially supported by Rus¬ 
sian Foundation for Basic Research. Yu.E.L. is supported 
by Program of Basic Research of HSE. 
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